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Abstract. The analogy of the nonlinear dislocation theory in crystals and the elec- 
tromagnetism theory is studied. The nature of some quantities is discussed. 



1. Burgers vectors and the Burgers space. 

Dislocations are one-dimensional defects of a crystalline grid used to explain 
the plasticity in crystals. Each dislocation line is characterized by its Burgers 
vector b, while the dislocated medium in whole is characterized by the so-called 
incompatible distorsion tensor T (see |TJ) . The components of the Burgers vector 




Fig. 1.1 Fig. 1.2 

for a dislocation line are determined by the following path integral along some 
closed contour encircling this dislocation line (see Fig. 1.1): 



(1.1) 



The Burgers vector b with components (1.1)| is teated as a vector of a special 
space, it is called the Burgers space. The Burgers space is an imaginary space, it 
is assumed to be filled with the infinite ideal (non-distorted) crystalline grid. The 



tensor T in (1.1) is a double space tensor: its upper index i is associated with some 
Cartesian coordinate system in the Burgers space, its lower index j is a traditional 
tensorial index associated with some coordinates y 1 
or curvilinear) in the real space where the crystalline medium moves 



y 2 , y 3 (no matter Cartesian 
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In the continual limit, when the number of dislocation lines is macroscopically 
essential, separate dislocation lines are replaced by their distribution p (see [l] for 
more details). Then (1.1) is replaced by the following integral equality: 



as j=1 



3 



pi v? ds. 



(1.2) 



Here S is some imaginary surface within the medium, n? are the components of the 
un it nor mal vector to S, and 7 = dS is the boundary of S (see Fig. 1.2). According 
to (1.2) , the double space tensorial quantity p is in terpreted as the Burgers vector 



per unit area. Applying the Stokes formula to (1.2) , we get the differential equality 

p = rotf. (1.3) 
Apart from (1.3) , we have the following equality (see [l]) : 



j = — grad w — 



dT 

~at' 



(1.4) 



The double space tensorial quantity j is interpreted as the Burgers vector crossing 



the unit length of a contour per unit time due to the moving dislocations (see [1 



However, the interpretation of the quantity w was not clarified in [1], This is the 
goal of the present paper. For this purpose below we study two special cases. 



2. Plastic relaxation. 

Let's consider a two-dimensional model of a crystalline medium with square 
cells (see Fig. 2.1). On the preliminary stage the crystal was distorted as shown on 



the Burgers 



space 

























































the real space 



Fig. 2.1 





Fig. 2.3 



Fig. 2.2. This distorsion is described by the following deformation map: 

x 1 =x 1 (y 1 ,y 2 )=y 1 -y 2 , 
x- = x z (y\y z ) = y\ 

Here we assume that x , x 2 are Cartesian coordinates in the Burgers space and 
y , y 2 are Cartesian coordinates in the real space, both are associated with the 
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orthonormal bases ei, e2 and Ei, E2 respectively. By differentiating (2.1) we find 
the components of the compatible distortion tensor T (see [if ): 



k dy k 



1 -1 
1 



(2.2) 



They are constants since the further evolution of our crystal goes without the 
displacement of atoms (see Fig. 2.3, Fig. 2.4, Fig. 2.5, Fig. 2.6). Hence, we have 



(2.3) 



Initially, our crystal has no dislocations at all. Therefore, the compatible and 
incompatible distorsion tensors are initially equal to each other: 



(2.4) 



1 t=t 



On Fig. 2.3 three pairs of the edge dislocations arise. Their Burgers vectors are 
b = ei and b = — ei in each pair. Therefore, the total Burgers vector of this group 
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Fig. 2.4 



Fig. 2.5 



Fig. 2.6 



of dislocations is equal to zero. On Fig. 2.4, Fig. 2.5, Fig. 2.6 the dislocations with 
negative Burgers vectors b = —ex move to the right. During the evolution time the 
blue arrow of the length 3 (see Fig. 2.3 and Fig. 2.6) is crossed by three dislocations 
with total Burgers vector b = —3 ex. The green arrow of the same length is not 
crossed by the moving dislocations at all. Therefore, we have the equality 



Jk 



dt 



t 



-1 




(2.5) 



at the center of our crystal. Behind the moving dislocations we find the undistortcd 
cells, they are marked by yellow spots on Fig. 2.6. This means that 



Tl 



1 
1 



(2.6) 



at the center of our crystal. Combining [Z2| |(23jt |(2~4)| , |(2~5)1 , and |(2~6)| we get 



j dt = 0. 



(2.7) 
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In the continuous limit, when the moving dislocations form the homogeneous and 
constant flow the above equality (2.7) can be transformed to the following one: 



dT . n 



(2. 



Comparing (2.8) with (1.4), we conclude that v = and T = const implies w = 
in our first example. 



3. Frozen dislocations. 

As the second example, we consider a three-dimensional crystal where the dis- 
location lines move together with the medium like water-plants frozen into the ice. 



the Burgers 
space 




v^'m/najy evolution map 

the interspace map 




Fig. 3.1 

The choice of Cartesian coordinates x 1 , x 2 , x 3 in the Burgers space is obligatory 
(see |jTj|) . In the real space we could choose either Cartesian or curvilinear coordi- 
nates. Below we choose Cartesian coordinates y , y 2 , y 3 for the sake of simplicity. 
Then the interspace map and its inverse map are given by the following formulas: 



f y 1 


^x 1 






y 2 


= x 2 


, y 3 


= X 3 



y 2 , 
y 3 - 



(3.1) 



The evolution of the crystal is subdivided into two stages (see Fig. 3.1). In the first 
stage, which is a preliminary one, the dislocations are produced: 



/ — 1 — 1 / ~1 ~9 

y = y (T,y ,y ,y a ), 
f = y 2 {-r,y l ,y 2 ,y 3 ), 
L f = y 3 (T,y\y 2 ,y 3 ), 



y 1 = y 1 ^^ 1 ,^ 2 ,^ 3 ), 
y 2 = ^(-r.y 1 ,^ 2 ,^ 3 ), 
y 3 = fi^y 1 ^ 2 ^ 3 ). 



(3.2) 
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In the second stage the dislocations are frozen and move together with the medium: 

(3.3) 



' V 1 = V 1 (t,T,y 1 ,y 2 ,y 3 ), 
y 2 = y 2 0,i-,y\y 2 ,y 3 ), 

. V 3 = 2/ 3 (*,T,y\£/ 2 ,|/ 3 ), 



y 1 = y 1 (t,T,y 1 ,y 2 ,y 3 ), 
y 3 = y 3 (i,T,2/\y 2 ,2/ 3 ). 



The compatible distortion tensor T at the time instant t is determined by the 
composite map including all of the three above maps |(3.1) , (3.2), and |(3.3) : 



(see [2J2] for comparison). Similarly, at the time instant t = r we have 



dx % 



(3.4) 



(3.5) 



From (3.4) and (3.5] , applying the chain rule to (3.3), we derive the relationship 

(3.6) 



9=1 



The relationship (3.6) expresses the evolution rule for the compatible distortion 
tensor T. If the dislocations are frozen into the material, then the incompatible 
distortion tensor T should obey the same evolution rule: 



n(t,y\y 2 ,y 3 



Y,T*(T,y\f,f) 

9=1 



dy^ 
dy k 



For the sake of simplicity, in the further calculations we denote 



t="3/'j- i 2 3-\ ®y q 

l&t,T,y ,y ,y) = Qyk' 



(3.7) 



(3.8) 



The quantities (3.8) form a non-degenerate square matrix T. Let S = T _1 be 

the inverse matrix for T and let Si = Sl(t, r, y 1 , y 2 , y 3 ) be the components of this 

inverse matrix. Note that T% are not the components of a tensor field. They arc 

not the components of a double space tensor in the sense of [1] as well. We shall 

fq 



not discuss the tensorial properties of the quantities T, 9 , we shall use them only as 



the notations for the partial derivatives (3.8) 



The following functional identity with two parameters t and r is quite obvious 
for the pair of mutually inverse maps given by the functions (3. 3)| : 



y k {t,T,y 1 {t,T,y 1 , ... ,y 3 ), ... ,y 3 (t,T,y\ ... ,y 3 )=y k . 



(3.9) 



By differentiating (3.9) with respect to the time variable t we easily derive 



r— 1 



(3.10) 
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Here v 1 , v 2 , v 3 are the components of the velocity vector v = vjt, y 1 , y 2 , y 3 ) of a 
point of the medium. Applying the partial derivative d/dy k to (3.10) and taking 
into account (3.8)| , we derive another useful equality: 



r=l " 

As for the equalities |(3.6) and (3.7), now they are written as follows: 
n = j2T^y\f,f)fl 

9=1 

n = jZK(r,y\y\f)n. (3.12) 

9=1 

Theorem 3.1. In the case of frozen dislocations the density of the Burgers vector 
p obeys the following evolution rule: 



p\ = dct f ]T £ 4( T > f) 9 qr S p r g pk . 

p—1 q—1 r—1 



(3.13) 



Here g p k and g qr are the components of the metric tensor and the dual metric tensor 
respectively (see [2] for details). 

Proof. The proof is pure calculations. The density of the Burgers vector is defined 
by the formula |( 1.3) . In coordinate form this formula is written as 



3 3 3 



r=l p=l q=l y 



(3.14) 



where uj rpq are the components of the so-called volume tensor (see [2]). The quan- 
tities <jf vq in (3.14] and the quantities g qr and g p k in |(3. 13) are constants because 
we chose the Cartesian coordinates y 1 , y 2 , y 3 (see Fig. 3.1). From |(3~. 12) we derive 



dT 



p ^ ^ 



dT* m {T,y\y 2 ,y 3 ) 
dyP 



dT? 



(3.15) 



Then we apply the chain rule to the first term and the formula (3.8) to the second 
term in the right hand side of the equality (3.15) . As a result we get 



3 3 



f)T l 

QyP Zs Is QyU 

a rn=l n=l " 

Applying the formula |(3.8) again, we derive 



dT l m {T,y\y 2 ,f) dy n 



-^-r + VPfr y 1 y 2 y 3 ) 

Q y p ^9 T "-rnKhy , </ ,y ) Q y p Q y q- 



~dyP ~ 2-^ 2s 

a m=l n=l 



dT; n ( T ,y\y 2 ,y 3 

Qy-n 



d 2 y m 
dyP dy q 
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Now we substitute the above equality into (3.14). The second term in its right hand 
side is symmetric in p and q. It vanishes when substituted into (3.14) because of 
the skew symmetry of uj rpq (see p| ). For p l k now we get 



3 3 3 3 3 cj'T i 2 / —1 —2 — 3\ 

A \\\\\^ , ,rpq al m\ T ->V iVjV) rfm <fr> 

Pk = l^2^l^Z^1^3kr^ pq ^ T p T q 



(3.16) 



r— 1 p—1 q—1 m—1 n—1 

In the next step we use the following well-known identity: 

3 3 3 



detT- ■ Jnm 



} rpq rpl rpn rpra 

r—1 p—1 q—1 

Since S = T^ 1 , from this identity we immediately derive 

3 3 3 



det TJ2 S 1 U lnm = E E ^ T p T 9 

p=l q=l 



(3.17) 



1=1 



Comparing (3.16) and (3.17) we can write the following formula for p k : 

3 3 3 3 rv ™ ' -,i -2 -3\ 



Pk 



r—1 I— 1 m—1 n—1 



r lnm 9T^(t , y 1 , y 2 , y A ) 

Lu 1 , 

dy n 



Now let's rewrite the formula (3.14)| for the time instant t = r: 

.mm dfi{r,y\f,f) 



3 3 3 

Pq(r,f,f,f) =EE E 9ql 



Z=l n—1 m—1 



Comparing (3.18) and (3.19) we can rewrite (3.18) as follows: 
p\ = det T £ E E ^ ST s' 9 (r, y\f,f)- 

r=l 1=1 9=1 



(3.18) 



(3.19) 



(3.20) 



Now it is easy to see that, in essential, |(3.20) coincides with the equality (3.13) 
which we needed to prove. So the proof is over. □ 

In the case of frozen dislocations the motion of the dislocation lines is com- 
pletely determined by the motion of the medium. Therefore, j should be expressed 
through v and p. In order to find this expression let's remember that j by definition 

is the total Burgers vector of the moving 
dislocations that cross the unit length of a 
contour 7 per unit time. It is clear that 
all of the dislocations passing through the 
dark parallelogram on Fig. 3.2 will cross 
the segment t during the next time inter- 
val dt. The total Burgers vector of such 
dislocations is determined by formula 




Fig. 3.2 



db l = ^p\n k dS, 



(3.21) 



fc=i 



s 
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where ndS = [v, r] dt and [v, r] is the cross product (vector product) of v and r. 
On the other hand, the same Burgers vector is given by another formula 



(3.22) 



c=l 

If we write the equality n dS — [v, r] dt in coordinate form 

3 3 3 



n k dS = J2J2J2 9 ks "sic v 1 t c dt, 



s=l 1 = 1 c=l 



then from |(3.21)| and |(3.22)| we derive the following formula for j : 

3 3 3 



j; 



(3.23) 



s=l 1=1 k=l 



Theorem 3.2. In th e case of frozen dislocation p and'} are related to each other 
by the formula |(3.23) . 



In order to continue the above calculations, now we substitute (3.16) into (3.23) 
As a result for the components of j we derive the following expression: 



3 3 3 3 3 3 
s=l 1=1 p=l q=l m— 1 n— 1 

Applying the well-known identity 1 



j , ,s Pq fl%»(T»y 1 .i/ a » ir > ) ^„ 



or 



nrm rprn 
P q 



to this expression we can decrease the multiplicity of summation symbols in it: 

dfl n {r,y\f,f)_ 

(3.24) 



3 3 3 

j^EEE^""^ 

p—1 m—1 n— 1 

-EEE^> ?f ; 

q—1 m—1 n—1 



dy n 



dr 



Theorem 3.3. In the case of frozen dislocation] and T are related to each other 
by the formula |(3.24) . 



Now we are going to calculate the time deriva tive of T using (3.12) . Upon doing 
this, we will be able to verify the equality (1.4) and calculate w: 



dt 2-*> 

q=l 



dT*(T,y\f,y^ 



dt 



f? + Y,r q (T,y\f,f)?§- 

q=l 



1 Here <5^, 51, 5^, and 5* are the Kronecker symbols. 
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g— 1 r— 1 

3 3 3 Qrfii/ -1 -2 -3\ 



9=1 P=l 

3 3 



EEE g 1: - p r ; T * -EE ^ 2 > 



In the above calculations we used t he ide ntities |(3.10)| and |(3.11)| . In order to make 
the resulting expression similar to (3.24) we change some summation indices: 



3 3 3 



p—l m—1 n— 1 



3 3 



P 

3 3 



E E^ 1 ,^ 3 ) £ ^ - E t^,v\M - q c 

m—1 q—1 m=l g=X 



dy c 



(3.25) 



Now let's return to the equation (1.4) . In coordinate form it looks like 



dTl 



dw l 



+ 7 = - - 

dt Jc dy c 



(3.26) 



Substituting (3.24) and (3.25) into (3.26), we derive 



dw - = E E E^" v9 fm ^^ T,yl,y2,y3s> 

q—1 m—1 n—1 



dy 



dy r > 



3 3 



m=lg=l y m =lg=l 



(3.27) 



Let's remember the formula (3.8) and apply it to T™ in (3.27) 



— = Vrfv dp m (r,y\y 2 ,y d ) dy 

QyC 2^ 2^ 

a m=\q=\ \n=l 



9y 



1 ,"2 ,-3^ fl,-n\ 



3 3 o g 3 3 f\rjrm 

E E S 8 ) |r^ m + EE ^(-, v\v\ f) v« 

m—1 q—1 m—1 q—1 



dy c 



(3.28) 



Looking at the right hand side of (3.28) , o ne eas ily recognizes the partial derivative 
of the product of three terms. Therefore, (3.28) is written as 



3 3 



£ = |;(EE^(^^ 3 )^ 
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Now, if we remember the formula (3.12), we obtain an even more simple equality: 

3 \ 



dw l 
dif 




(3.29) 



Integrating the equality (3.29) , we define the components of the vector w: 



9=1 



(3.30) 



This definition is unique up to some inessential terms depending only on the time 
variable: w' — > w l + w l (t). These terms can be omitted since w is used only in the 
form of its gradient (see the equation (1.4)). 



Theorem 3.4. In the case of frozen dislocation the vectorial parameter w is de- 



termined by the formula (3.30) 



(3.13) 


, (3.24) 


ions |(1.3) 


and 



4. The analogy to electromagnetism 
and the gauge transformations. 



The basic differential equations (1.3) and (1.4) describing the kinematics of dislo- 
cations in crystals are similar to the equations expressing the electric and magnetic 
fields E and H through the scalar potential <p and the vector potential A: 



H = rot A, 



E 



grad (p 



1 dA 
c ~dt' 



(4.1) 



These differential equations (4.1) admit the following gauge transformations 

1 dtp 



A — > A + grad-0, 



<P Q , 

c at 



(4.2) 



that change potentials, but do not change the actual physical fields E and H (see 
p] for the reference). By analogy we can write the gauge transformations 



T + gradV, 



dtp 
~dt 



(4.3) 



for T and w in th e di fferen tial equations (1.3) and (1.4) . Despite to the striking 
similarity of (4.2) and (4.3)| their roles arc absolutely different. The matter is that 
A and ip in the electromagnetism are not actual physical fields, they are derived 
mathematically from the Maxwell equations (see |3|). Unlike A, the incompatible 
distorsion field T is an actual physical field describing the state of a medium. 
Applying the gauge transformation (4.3)| to it, we get another field describing some 
other state of the medium. For this reason, the nonlinear theory of dislocations is 
not a gaug e theo ry or, at least, it is a gauge theory with the different gauge group, 
other than |(4.3) . 



GAUGE OR NOT GAUGE ? 



11 



The physical nature of the parameter w is still misty. The following conjecture 
opens a way to clarify it. It is based on the theorem 3.4. 



Conjecture 4.1. The vectorial parameter w in the equation (1.4) is always deter- 
mined by the formula |(3.30) . 



If we admit the conjecture 4.1, then the differential equation (1.4) describing the 
time evolution of the incompatible distortion tensor T can be written as 



&r, 



(4.4) 



where V p = d/dy p since we use the Cartesian coordinates y 1 , y 2 , 



and 



rppi 



(4.5) 



In this form (4.4) the differential equation (1.4) is quite similar to the balance 
equations which are traditional in the mechanics of continuous media (see the mass 
balance equation, the momentum balanc e equ ation, and the energy balance equa- 
tion in [4] for comparison). The tensor (4.5) is interpreted as the density of the 



incompatible distorsion flow. 

A conjecture is not a theorem yet. 
possible options for the parameter w in the equation (1.4) 



Therefore, we need to discuss all of the 



(1) the parameter w is determined by the formula (3.30) 

(2) the parameter w is determined by the formula 



dt 



(4.6) 



(3) the parameter w is determined by some formula other than (3.30) and (4.6' ; 

(4) the parameter w is an independent parameter of a dislocated medium. 

The option (1) is my favorite option. It is supported by the above two examples 
considered in the sections 2 and 3. 



The option (2) is often chosen in many papers. Substituting (4.6) into the 



equation 



1.4), we obtain the following equality: 



dT _ dT 
~dt ~~dt 



(4.7) 



The equality (4.7)| is referred to as the additive decomposition of the total distorsion 
into the elastic and plastic parts. It is consistent in the linear theory, where both T 
and T are approximately equal to the unit matrix. In the nonlinear case the equality 
1(4.6) is not a valid option since it contradicts the theorem 3.4. The multiplicative 
decomposition of the total deformation tensor suggested in [[4] is more preferable. 

The choice of the option (3) or the option (4) produces more questions than the 
answers. One should carry out a special research in order to exclude both these 
options or to prove one of them. 
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